Multilevel Tunnelling Systems and Fractal Clustering 
in the Low- Temperature Mixed Alkali-Silicate Glasses 
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The thermal and dielectric anomalies of window-type glasses at low temperatures (T < 1 K) 
are rather successfully explained by the two-level systems (2LS) standard tunneling model (STM). 
However, the magnetic effects discovered in the multisilicate glasses in recent times, magnetic effects 
in the organic glasses and also some older data from mixed (Si02)i-a:(K20)a; and (Si02)i-i!(Na20)a: 
glasses indicate the need for a suitable extension of the 2LS-STM. We show that - not only for the 
magnetic effects, but already for the mixed glasses in the absence of a field - the right extension of the 
2LS STM is provided by the (anomalous) multilevel tunnelling systems (A-TS) proposed by one of us 
for multicomponent amorphous solids. Though a secondary type of TS, different from the standard 
2LS, was invoked long ago already, we clarify their physical origin and mathematical description 
and show that their contribution considerably improves the agreement with the experimental data. 

PACS numbers: 61.43.Fs, 77.22.-d, 77.22.Ch, 65.60.-fa 
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I. INTRODUCTION 



Glasses are ubiquitous materials of considerable im- 
portance for many practical applications, however for 
physicists the nature of the glass transition and the ulti- 
mate microscopic structure of glasses determining their 
physical properties remain to this day issues of consid- 
erable intellectual challenge 0- Glasses are normally 
regarded as fully homogeneously disordered amorphous 
systems, much alike liquids except for the glassy arrested 
dynamics close and below the glass transition tempera- 
ture Tg, which leads to an increase of several orders of 



magnitude in the viscosity for T — >■ T+. Nevertheless 
this homogeneity is most probably only a useful ideal- 
ization, for real glasses must always contain some small 
(in ceramic glasses not so small) concentration of tiny, 
ordered or nearly-ordered regions of variable size with 
their own frozen dynamics. Indeed the thermodynam- 
ically stable phase of an undercooled liquid would be 
the perfect crystal, thus every substance in approaching 
the crystallization temperature Tc (Tc > Tg) from above 
would spontaneously generate local regions of enhanced 
regularity (RER) much like a system (a vapour or a para- 
magnet) approaching its critical temperature is known 
to develop regions (droplets) resembling the ordered low- 
temperature phase. These RER are of course to be dis- 
tinguished from the concept of short-ranged atomic order 
which is typical of ideal glasses and is restricted to the 
first few atomic spacings. We are considering in this pa- 
per realistic glasses in which a degree of devitrification 
has occurred. The size and concentration of these RER 
will depend, e.g., on the rapidity of the quench leading 
to the formation of the glass, but also on the chemical 
composition of the substance, the presence of impurities 
and so on. However, on general grounds, even the purest 
of glasses should contain RER in non-zero concentration 



and size. 

That this is the case has been demonstrated recently 
for the structure of the metallic glass Zt^qCvl^^AI^ [3], 
where a combination of fiuctuation electron spectroscopy 
(FEM) and Monte Carlo simulation (MC) has revealed 
the presence of crystalline regions of sub-nanometer size 
embedded in an otherwise homogeneously amorphous 
mass of the same composition. It is believed that other 
metallic glasses should present similar structural fea- 
tures and thus - on general grounds ~ one would ex- 
pect that non-metallic window glasses too, like pure Si02 
and all the more so the commercial multisilicates of com- 
plex chemical composition, should present a multiphased 
structure with the size and concentration of the near- 
crystalline regions, or RER, depending (e.g.) on compo- 
sition, quench rate and the presence of impurities act- 
ing as nucleation centres for the RER. Indeed, materials 
of the general composition (MgO)2;(Al2 03)y (Si02)i-x-y 
(MAS, in short) are termed ceramic glasses (one of the 
best-known commercial examples being Schott's Ceran 
where Li2 replaces MgO, and of course CaO or BaO 
can also replace or be added to MgO and still yield 
a ceramic glass). These materials are known to con- 
tain micro-crystals embedded in an otherwise homoge- 
neously amorphous matrix Q. This is not surprising, 
for materials made up of a good glass-former (e.g. Si02, 
AI2O3 ...) and good crystal-formers (e.g. BaO, K2O ...) 
are known to be multiphased [5i] with the good crystal- 
formers generating their own pockets and channels carved 
out within the otherwise homogeneously amorphous net- 
work of the good glass-former's ions J|. Within these 
pockets and channels, incipient nano- or even micro- 
crystals may form, but the point of view will be taken 
in this work that on general grounds even the purest, 
single-component (e.g. As, Si02) glass- former will be 
rich in RER unless the quench-rate from the melt is so 
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large as to avoid the formation of crystalline regions or 
RER. 

These refined structural details of glasses are evidently 
hard to reveal in all and especially the near-ideal cases 
(no good crystal-formers, no impurities added, and rapid 
quenches) with the available spectroscopic techniques. 
For example X-ray spectroscopy does not reveal nano- 
crystals below the nanometer size. However, at low- and 
very low-temperatures - where all said structural fea- 
tures remain basically unaltered - some recent exper- 
imental findings might now improve perspectives with 
what would appear set to become a new spectroscopy 
tool. Indeed a series of remarkable magnetic effects have 
recently been discovered in non- mag netic glasses (mul- 
tisihcates and organic glasses) [7-ll| with, in the opin- 
ion of the present Authors, a most likely explanation for 
the new phenomena stemming precisely from the multi- 
phase nature of real glasses and the presence of the RER 
or micro-crystalline regions in their microscopic struc- 
ture. In turn, when the multiphase theory shall be fully 
developed, the magnetic effects could represent a valid 
new spectroscopic tool capable of characterizing micro- 
or nano-crystals or even incipient crystals and RER in the 
real glasses. The key to this possible development is some 
new exciting physics of the cold glasses in the presence 
(and even in the absence, as shown in the present pa- 
per) of a magnetic field. The magnetic effects in the cold 
glasses could become, eventually, the amorphous coun- 
terpart of the de Haas - van Alphen and Shubnikov - de 
Haas effects in crystalline solids in determining the real 
structure of amorphous solids. 

Systematic research in the low-temperature properties 
of glasses has been on-going for more than 40 years and 
some significant theoretical and experimental progress 
has been made in the understanding of the unusual be- 
haviour of glasses and of their low-temperature anomalies 
[l^ . This temperature range (T <1 K) is deemed impor- 
tant for the appearance of universal behaviour (indepen- 
dent of composition), as well as for the effects of quantum 
mechanics in the physics of glasses. However, to make 
progress in the understanding of the low-temperature 
physics of glasses there remains a wide range of important 
questions that are still open or only partially answered, 
particularly in the light of some still poorly understood 
recent, and even older, experiments in cold composite 
glasses. 

It is well known that cold glasses show somewhat uni- 
versal thermal, acoustic and dielectric properties which 
are very different from those of crystalline solids at low 
temperatures [l^. Below 1 K, the heat capacity Cp of 
dielectric glasses is much larger and the thermal conduc- 
tivity K is orders of magnitude lower than the correspond- 
ing values found in their crystalline counterparts. Cp de- 
pends approximately linearly and k almost quadratically 
on temperature T, whereas in crystals one can observe 
a cubic dependence for both properties, well understood 



in terms of the Debye-Einstein's theory of lattice vibra- 
tions. The dielectric constant (real part) e' and sound 
velocity at low frequencies display in glasses a univer- 
sal logarithmic dependence in T. These "anomalous" 
and yet universal thermal, dielectric and acoustic prop- 
erties of glasses are well explained (at least for T < 1 K) 
since 1972 when Phillips [1% and also Anderson, Halperin 
and Varma [15j , independently, introduced the tunnelling 
model (TM), the fundamental postulate of which was 
the general existence of atoms or small groups of atoms 
in cold amorphous solids which can tunnel like a sin- 
gle quantum-mechanical particle between two configura- 
tions of very similar energy (two- level systems, 2LS). The 
2LS TM is widely used in the investigation of the low- 
temperature properties of glasses, mostly because of its 
technical simplicity. In fact, it will be argued in this 
paper that tunneling takes place in more complicated 
local potential scenarios (multiwelled potentials) and a 
situation will be discussed where the use of a number 
of "states" greater than two is essential. Moreover, new 
insight will be given on the role of percolation and frac- 
tal theory in the TM of multicomponent glasses. We 
present in this paper the justification and details of the 
construction of an extended TM that has been success- 
fully employed to explain the unusual properties of the 
cold glasses in a magnetic field [3], as well as in zero 
field when systematic changes in the glass' composition 
are involved 



17 . 



The linear dependence in ln(T) of the real-part of the 
dielectric constant e'{T) makes the cold glasses useful in 
low-temperature thermometry and, normally, structural 
window-type glasses are expected to be isotropic insula- 
tors that do not present any remarkable magnetic-field 
response phenomena (other than a weak response in Cp 
to the trace paramagnetic impurities). For some multi- 
component silicate glass it has become possible to mea- 
sure observable, much larger than expected changes in 
e'{T,B) {Se'/e' ^ 10~^) already in a magnetic field as 
weak as a few Oe Q. A typical glass giving such strong 
response has the composition Al203-BaO-Si02: thus a 
MAS ceramic-glass, herewith termed AlBaSiO. The mea- 
surements were made on thick sol-gel fabricated films, 
a fabrication procedure favoring micro-crystal formation 
0, cooled in a '^He-'^He dilution refrigerator reaching 
temperatures as low as 6 mK. Magnetic effects have been 
reported for both the real- and imaginary-part of e at 
low frequency (w ~ 1 kHz), for the heat capacity Cp (see 
e.g. [16|) and for the polarization echo (where changes in 
the presence of a magnetic field have been the strongest 
[lo| ) as well. This behavior was confirmed in other mul- 
ticomponent glasses, like borosilicate optical glass BK7 
and commercial Duran ^ and, moreover, similar effects 
on e'(T) have been confirmed in studies of the struc- 
tural glass a-Si02+xCj;H^ in the range 50< T <400mK 
and B <3T Although the dielectric magnetocapaci- 
tance enhancement is not dramatic (6e'{B)/e' is typically 
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in the 10~^ - lO"** range), the available measurements 
show that an unusual effect of the magnetic field is in- 
deed present in the above glasses, yet not measurable in 
ultra-pure Si02 (Suprasil W), and cannot be ascribed to 
spurious agents 18[ or to trace paramagnetic impurities 
(always present in silicate glasses, although in <6 ppm 
concentration in the case of BK7). Polarization-echo ex- 
periments in the AlBaSiO, Duran and BK7 glasses have 
also shown considerable sensitivity in the response of the 
echo amplitude to very weak magnetic fields and the mag- 
netic effects clearly do not scale with the concentration 
of paramagnetic impurities [1, [13] ■ Striking magnetic ef- 
fects, the presence of a novel isotope effect and remark- 
able oscillations in the dephasing time have also been 
reported in studies of the polarization echos in organic 
glasses (amorphous glycerol) \v\\. However, in terms of 
a detailed theoretical justification for all of the observed 
magnetic effects (and the lack of an observable magnetic 
effect in the acoustic response so far) an explana- 
tion relying on a single theoretical model for all of the 
available experimental data is still missing. We believe 
the two-phase model reproposed in this paper to be the 
correct generalization of the standard 2LS TM that is be- 
ing sought and here we work out its predictions in zero 
magnetic field, but for different controlled concentrations 
of glass-forming and crystal-forming components. In this 
way, we put our approach to a new test. 

The essential behavior of the dielectric response of 
glasses at low temperatures is well known [13|: start- 
ing from the lowest temperatures, the dielectric constant 
e' first decreases with increasing temperature due to the 
resonant interaction of the electric field with the tunnel- 
ing systems (TS). According to the standard 2LS TM 
(STM from now on), the dielectric constant is predicted 
to vary like -In T due to the constant density of states of 
the TS. Above a certain temperature To{uj) relaxational 
absorption of the TS becomes important, resulting in an 
increase of the dielectric constant with temperature pro- 
portional to +lnT according to the STM. This has been 
checked experimentally for several glasses. The temper- 
ature Tq of the resulting minimum depends on the fre- 
quency oj and occurs around 50 to 100 mK in measure- 
ments at about 1 kHz: in AlBaSiO glass, the position of 
the minimum has been shown to shift considerably when 
the amplitude of the driving electric field is considerably 
increased (probably, the consequence of non-linearity set- 
ting in). 

Some more interesting behavior has been shown by 
some as yet unexplained data from experiments on the 
mixed (Si02)i-a;(K20)2: and (Si02)i-2:(Na20)2; glasses, 
studied as a function of the concentration x of the good 
crystal- former at low temperatures [2^ . The heat capac- 
ity Cp{T) for these glasses is larger than that for pure 
vitreous silica and the behavior as a function of T is very 
peculiar for different molar concentrations x of potassium 
or sodium oxide and is not explained by the STM. The 



heat capacity decreases and then increases again with in- 
creasing molar concentration x of K2O. The minimum in 
the dielectric constant e'(T) is observed for Tq near 100 
mK as is typical for these glassy solids. The tempera- 
ture dependence of e', both above and below Tq, shows 
however a slope in ± In T qualitatively increasing with in- 
creasing concentration x of K2O. One can notice, more- 
over, that above the minimum Tq the relaxation part 
of e' is increasing faster in slope than the resonant part 
below To for the same x [20| . a feature completely unex- 
plained thus far. This work is an indication that not only 
the magnetic and electric fields influence the properties 
of glasses, but the concentration of chemical species in 
the composite materials too (a fact not accounted for by 
the STM). In this paper we show in detail how the very 
same approach that explains the magnetic properties in 
the multisilicates [3] also provides a quantitative expla- 
nation for the above-mentioned composition-dependent 
physical properties. The picture that emerges regarding 
the nature of the TS in the multicomponent glasses pro- 
vides a novel and detailed description of the micro- and 
nano-structure of the glassy state. In turn, the linear 
dependence of the concentration xats of anomalous TS 
(ATS) - those responsible for the magnetic and compo- 
sition effects in our theory - on a; fully corroborates the 
founding assumptions of our approach. 

The paper is organised as follows. In Section II we 
present a detailed justification for the two-phase ap- 
proach and the construction of the two-species TS model 
for the amorphous solids at low temperatures. In Section 
HI we present the detailed predictions of this model for 
the dielectric constant e'(T, x) as a function of temper- 
ature T and composition x of alkali oxide (good-crystal 
former) for the mixed glasses and we compare the predic- 
tions with the experimental data |20|]. In Section IV we 
present the detailed predictions of our model for the heat 
capacity Cp(T, x) for the mixed glasses and we com par e 
the predictions with the available experimental data I20j . 
Section V contains our conclusions about the nature of 
the TS, namely we show how the tunneling "particle" 
must in fact represent a whole cluster of N correlated 
real tunneling ions in the material. Finally in the Ap- 
pendix we work out how the effective tunneling param- 
eters of our model are related, via iV, to more standard 
microscopic tunneling parameters. A short preliminary 
account of this work was published in [17 |. 



II. BUILDING UP A SUITABLE TUNNELING 
MODEL 



The traditional picture [13| viewed the TS, present in 
low concentration ('^ 10"'^^ g~^) in the material, associ- 
ated with the non-equivalence of two (or more) bonding- 
angle configurations per atomic unit in the amorphous 
solid's atomic structure. Each TS is represented in 
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the standard case by a particle in an asymmetric (one- 
dimensional (ID)) double-well potential where, at low 
T, only the ground states of the two constituent single 
wells are assumed to be relevant. Consequently, only 
the two lowest-lying double-well states are taken to de- 
termine the physics of each single TS. A 2LS simplified 
picture then applies and one can describe the low en- 
ergy Hamiltonian of each independent TS in terms of 
an equivalent notation with spin-^ pseudospin matrices 
ax and Uz (Pauli matrices), leading to the compact no- 
tation h'^^ = — i(A(Tz -|- ApCTa;) for the Hamiltonian 
of a single 2LS TS. In matrix form (the so-called well- 

or position-space representation, (il-ffg N) being the 
two unequivalent wells, i = 1, 2 or i = L, i?) this then 
reads: 



(2) 



A Ao 

Ao -A 



(1) 



Here, the phenomenological parameters A and Ao 
(known as the energy asymmetry and (twice) the tun- 
nelling matrix element, respectively) represent a way of 
describing the essential low-T relevant features of the full, 
and yet unknown in its details, TS single-particle Hamil- 
tonian in the effective single-well matrix representation. 
One obtains £i_2 — i^-^/A^ -l- Aq for the two lowest- 
lying energy levels and the physics of the glass is then 
extracted by assuming (initially) the 2LS to be indepen- 
dent entities in the glass and averaging physical quanti- 
ties over a probability distribution for the parameters A, 
Aq of the standard form (P being a material-dependent 
constant): 



P(A,Ao) = — . 



(2) 



This distribution reflects the generally accepted opinion 
that A and — ln(Ao/?iri) (the latter proportional to the 
double- well potential barrier Vq divided by the single- well 
attempt frequency VL, Vq/MI) should be rather broadly 
distributed in a homogeneously disordered solid. This 
leads to an almost constant density of states (DOS) and 
the above STM has been employed with considerable suc- 
cess in order to explain a wide range of physical proper- 
ties (thermal, dielectric (ac and pulsed), acoustic, and so 
on [13[) of non-metallic glasses below 1 K. 

There are, however, several drawbacks with the STM 
as thoughtfully pointed out by Leggett and co-workers 
21| . For a start, the nature of the TS (and of the two 
wells of a single 2LS) and that of the motion inside a sin- 
gle TS remain to date completely unknown [22|. Much 
easier is the diagnostic for the nature of 2LS in the case 
of disordered crystals, such as Li-KCl or KBr-KCN so- 
lutions [2^ (we shall come back to disordered crystals 
later). On general grounds, other types of (multilevel) 



excitations are always possible in glasses and it is not 
clear why their distribution of parameters should be so 
similar (and given by Eq. ([2])) in all of the amorphous 
solids. Next, the STM has gathered great consensus for 
the explanation of many experiments at low tempera- 
tures, but in its simplest form (Eqs. ([I])-©) it fails to 
explain sound velocity shift and adsorption data at low- 
T and the origin of the "bump" in Cp (and "plateau" in 
k) well above Tq that goes under the name of boson peak 
(see e.g. the references in [lH). Moreover, the STM fails 
to explain the remarkable universality of the ultrasonic 
attenuation coefficient (roughly, independent of ev- 
ery external parameter and glass chemical composition) 
below 1 K |26|. To resolve these (and other) difficulties 
with the STM, Leggett and collaborators have proposed 
a generic model in the context of anharmonic elasticity 
theory which can account for all of the significant features 
of glasses below 1 K, including the super universality of 

However, it is hard to see how this generic elastic 
model can be extended to account for the magnetic and 
composition-dependent effects in glasses, also considering 
that in the multicomponent (i.e. real, not model) glasses 
most of the said universality features (e.g. in Cv(T, B) 
and e'(T,B) 0, [l^ or in Cp(T,x) and e'(T,x) U^^) 
are lost. Therefore, here we adopt the strategy of resum- 
ing the TS approach by means of a completely different 
(and more modern) justification for the TM, and then 
extend the STM to take the presence of a magnetic field 
into account and to explain composition-dependent fea- 
tures (this work) . In a rather general fashion, the TS can 
be thought of as arising from the shape of the theoreti- 
cal potential-energy landscape E{{ri}) of a glass as T is 
lowered well below the glass freezing transition Tg. The 
concept of free-energy landscape was introduced, e.g., by 
Stillinger [g^] and successfully employed in the study of 
glasses (e.g. 0) and spin-glasses (e.g. [11]). A large 
number of local and global minima develop in E{{ri}) as 
T — )■ 0, the lowest-energy minima of interest being made 
up of riw = 2,3,... local wells separated by shallow en- 
ergy barriers. At low-T these configuration-space local 
multiwelled potentials are our TS and it seems reason- 
able to expect that the = 2 - welled potentials (2LS) 
should be ubiquitous in this picture. These should be 
thought of as an effective representation of local "trem- 
blements" of the equilibrium positions {r^''^} of some 
of the glass atoms/ions' positions spanning over a large 
number of near- neighbors' distances (unlike in the case 
of disordered crystals, where the TS are known to be 
rather well- localized dynamical entities). Hence, just as 
the = 2 - welled case is possible, so ought to be the 
Uw = 3, 4, ... - welled situations which would also be local 
rearrangements involving several atoms/ions/molecules. 
The concentration of these local potentials should not 
necessarily decrease exponentially with increasing , in 
glasses, as it is known to happen for the disordered crys- 
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tals (2LS present with probability c^, 3LS with c^, 4LS 
with ... and so on, c being the defects' percent con- 
centration) . 




FIG. 1: (color online) The energy landscape (for p —1 
Lennard- Jones density, adapted from Ref. 29|) of a toy 
glass model, with highlighted multiwelled potentials (black 
the 2LS, light blue the 3LS, 4LS, ...). 

We can reason this out over the quantitative descrip- 
tion of the glassy energy landscape of a model situ- 
ation, as was studied by Heuer [29| who considered 
the molecular-dynamics (MD) simulation data of a toy 
glass made up of several (13 or 32) particles interact- 
ing through a Lennard- Jones potential and with periodic 
boundary conditions applied. Adopting a suitable ID 
projection procedure, where a "distance" between two 
local total energy minima is (not completely unambigu- 
ously) defined, the ID position of a local minimum is 
somehow attained and the energy landscape of the model 
system can be charted out. Fig. [T] reports this chart for 
the total energy landscape for a given density (from [29|). 
Beside the deep minimum of the crystalline configuration, 
a large number of local minima is visualized and then a 
suitable definition of local double- welled potentials (2LS) 
is adopted to classify couples of adjacent minima consti- 
tuting a single tunneling 2LS (highlighted in black, in 
Fig. [T]). This definition guarantees that at low temper- 
atures a "particle" subjected to any such local poten- 
tials will switch between both minima without escaping 
to a third minimum. Interestingly, the distribution of the 
tunneling parameters A, Aq (suitably defined) for these 
2LS could also be evaluated from MD simulations of the 
above toy model, and this P(A, Aq) turned out to be 
not so perfectly flat as a function of A as implied by Eq. 
([2]). Rather, an increase (though no divergence) of prob- 
ability for 2LS with A — >■ was measured in previous 
MD simulations [s^ . Still, Fig. [T] also allows for tun- 
neling multiwelled local potentials to be identified, and 
we have highlighted (in light blue) some of them (three- 
and four- welled local potentials). The requirement that 
a "particle" subjected to such multiwelled local poten- 
tials should not escape (at low-T) to foreign minima has 
been equally respected and one can see that these multi- 
welled situations are not at all rare. We therefore believe 
that 3LS, 4LS and so on should also be considered in 



the TM. The reduced Hamiltonians (well- or position- 
representation) for these local multiwelled potentials can 
be easily written down, as generalizations of Eq. ((TJ. For 
rill, — 3 (3LS): 



^ \ Do E2 Do 

Do Do Es 



El Do Do 



(3) 



where Ei, E2, E^ are random energy asymmetries be- 
tween the wells chosen to satisfy Yl^i=i^i — ^^"^ taken 
from an appropriate probability distribution (see below) , 
together with the tunneling parameter Dq > (see be- 
low). For n^^A (4LS): 



H, 



(4) 



/ El Di L>2 Di \ 

Di E2 Di D2 

D2 Di E3 Di 

\ Di D2 Di Ei ) 



(4) 



where Ei,E2,E3,E4 arc random energy asymmetries 
taken from an appropriate probability distribution, to- 
gether with the tunneling parameters Di (n.n. well hop- 
ping) and D2 (n.n.n. hopping, \D2\ <^ \Di\)- These 
are simple, possible choices; clearly, other special-purpose 
generalizations of the 2LS matrix Hamiltonian are possi- 
ble and we believe that the 3LS of Eq. (|3]) is the minimal 
generic multiwelled potential wich can take the magnetic 
field into account [iq (the 2LS Hamiltonian of Eq. ([T]) 
could also be adjusted for this purpose, however the en- 
ergy spectrum would be totally insensitive to B). One 
can easily convince oneself, at this point, that so long as 
the energy parameters of the above multiwelled effective 
Hamiltonians obey the usual uniform distribution (Eq. 
dH), suitably reformulated) as is advocated by the STM, 
the DOS g{E) will remain (roughly) a constant. It is then 
to be expected that all these multiwelled local potentials 
will give rise to the very same physics as in the n^, — 2 
case and that thus, in practice, the 2LS choice repre- 
sents the appropriate minimal model for all of the extra 
low energy excitations characterising amorphous solids 
at low-T. It is clear from the above discussion, however, 
that the 2LS tunneling "particle" is no atomic particle 
at all, but, on general grounds, it rather represents the 
local rearrangements of a good number of real particles 
(ions or molecules). 

All changes if the glass is made up of a mixture of 
network-forming (NF) ions (like those of the good glass- 
forming Si04 or (A104)~ tetrahedral groups) as well as 
of network-modifying (NM) ions (like those of the good 
crystal-forming K+ or Na"*", or Ba^"*", ... from the rela- 
tive oxides) which, these last ones, could act as nucle- 
ating centres for a partial devitrification of the glass, 
as is known to occur in the multicomponent materials 
31 , l^l . Indeed the NM-ions of the good crystal- formers 



FIG. 2: (color online) Molecular dynamics snapshot of the 
structure of sodium trisilicate at 2100 K at the density p = 2.2 
g cm~^. The big blue spheres that are connected to each 
other represent the Na atoms. The Si-O network is drawn by 
yellow (Si) and red (O) spheres that are connected to each 
other by covalent bonds shown as sticks between Si and O 
spheres (from by permission). 



FIG. 3: (color online) A region including the crystal-like 
supercluster from a snapshot of the model simulation - in- 
corporating fluctuation electron microscopy data - of the 
Zr5oCu45Al5 metallic glass at 300 C (from Q]). The atomic 
separation distances of the middle zone are about 0.25 nm. 
This is a first realistic image of a crystal embryo in a glass; 
this object should not be confused with the concept of short- 
range order in ideal glasses. 



are termed "glass modifiers" in the glass chemistry liter- 
ature [33I since they do not become part of the intercon- 
nected random network but carve out their own pockets 
and channels within the glassy network [sj] ■ Fig- [H 
(courtesy from W. Kob, 6]) shows a snapshot of a MD 
simulation of the glass having composition Na20-3(Si02) 
(or, (Na20)o.25(Si02)o.75) at 2100 K (above Tg, in fact) 
in which the non-networking NM Na-atoms are put in 
evidence (big blue spheres). Simulations and experi- 
ments in the multisilicates definitely show that the NM- 
species in part destroy the networking capacity of the NF- 
ions and form their own clusters inside the NF-network 
0. The chance for these NM-clusters to be the nest of 
RER, incipient- or actual micro-crystals is obviously very 
good, considering that these clusters are made of good 
crystal-forming atoms. However, on general grounds 
and as discussed in the Introduction, we shall take the 
attitude that even the purest single-component glasses 
will contain RER in some measure. Fig. [3] (from Q) 
shows one such RER within a snapshot from a joint 
MC-simulation/FEM- measurement on the metallic glass 
Zr5oCu45Al5. The picture clearly shows an embryo crys- 
tal which could not grow to macroscopic size due to the 
arrested dynamics below Tg] such structures are expected 
to ubiquitous in all glasses, metallic and non-metallic 
35| . Except that they are difficult to observe with 
the available spectroscopic tools when sub-nanometric in 
size. The concentration and size of these RER will dictate 
whether magnetic- or composition-effects become mea- 
surable in the low-T experiments. a-Si02 in its purest 
form (Suprasil W) revealed no measurable magnetic ef- 
fects |lS[i3i. 

It goes without saying that TS forming in the prox- 
imity and within these RER or micro-crystalline regions 



will require a completely different mathematical descrip- 
tion, in particular the possibility of having more than two 
wells affords a more realistic description of the energy 
landscape. Hence, > 2 multiwelled systems inside 
the glass-modifying NM-pockets and -channels should 
follow some new energy-parameters' distribution form 
when some degree of devitrification occurs, leading to 
entirely new physics. One of the present Authors has 
proposed that precisely this situation occurs inside the 
magnetic-sensitive multicomponent glasses [l6| . and in 
this paper we show how this theory explains the B = Q 
composition-dependent dielectric and heat capacity data 
of [20| as well. Instead of the standard ID double- welled 
(W-shaped) potential, leading to Eq. ([T]), which contin- 
ues to describe the ordinary tunneling 2LS inherent to 
the homogeneously-disordered a-Si02 network, we take 
for the TS nested in or near the RER, crystal embryos or 
micro-crystals, the model of a "particle" having charge q 
and moving in a rim-welled 3D potential of the shape dis- 
played, for riw = 3, in Fig. |3]for the 2D (x, 2/)-space. The 
hopping Hamiltonian of a single, non-interacting tunnel- 
ing 3LS has therefore the form (for a fictitious second- 
quantization particle in the well-coordinate representa- 
tion) 



and is described in matrix form by Eq. ^ (where in fact 
{i\HQ is displayed, |i) (i =1,2,3) denoting the single- 
well ground states). This is our minimal generic model 
for a multiwelled TS. The parameter Dq is chosen posi- 
tive (contrary to custom in the STM, indeed — 5A0 < 
in Eq. ([1}) for a good number of reasons. First, due 
to the possible softness of the local NM-potential, since 
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indeed in general 13| Dq ~ aTifie"''^"/''^, a and b being 
numbers such that for Vq ^ hil a > and b = 0{1) can 
arise 13, [l^. This choice is still compatible with the 
concept of tunneling and at the same time yields rather 
large values of Dq « hQ. On more general grounds, how- 
ever, one should take into account that the tunneling 
"particle" is not moving in a vacuum, but is embedded 
in a solid that is for the most part deprived of micro- 
scopic dynamics, at low-T. Thus the surrounding frozen 
atoms are taking a part in the determination of the tun- 
neling particle's lowest stationary states. In the case of 
a perfectly CS-symmetric local n^^ — 3 welled potential 
of the type depicted in Fig. SI Hamiltonian ^ leads 
to a doubly-degenerate ground state and a first excited 
non-degenerate state (as is easily verified from Eq. ([3]) if 
El = E2 = E3). This may seem unphysical and yet Suss- 
mann has demonstrated, in a remarkable paper ^36,] , that 
for electrons trapped in a crystal (or equivalently in a 
glass) the situation above described is realised whenever 
the trapping potential is multiwelled with a triangular 
{uu, = 3) or tetrahedral (n^, = 4) well-centers geometry. 
The binding of the seemingly anti-bonding ground state 
is then guaranteed by the TS interaction with the rest of 
the solid. This reasoning is irrelevant for the STM-2LS 
parameter Aq, since both positive and negative signs for 
this parameter yield the same physics. If rti„ > 2 the sign 
will matter and Sussmann's work shows that the choice 
Do > is physically justified for an embedded parti- 
cle in the glass (or vacuum). Finally, it will be shown 
in the Conclusions that in fact the tunneling "particle" 
cannot be considered a single atom, ion or molecule, but 
rather it represents a cluster of N correlated tunneling 
atomic-scale particles, with N w 200. Then it is reason- 
able to expect that the ground state of such a cluster 
might be near-degenerate, so our choice -Dq > for the 
effective single tunneling "particle" is sound and not in 
conflict with any general quantum-mechanical principle. 
This Dq > is the major assumption for the multiwelled 
TS theory. 

At this point, we make a choice for the probability dis- 
tribution of the parameters Ei, E2, E3 and Dq of a tun- 
neling 3LS nesting in the proximity of a RER, crystal em- 
bryo or micro-crystal (one could also work with a = 4 
model potential, in the Appendix we show that essen- 
tially the same results can be attained). This is dictated 
by the fact that near-degeneracy {Ei = E2 = E3) must 
be favored, yet not fully attained for the wells' energy 
asymmetries of one such 3LS. We thus choose, assum- 
ing again the tunneling potential barriers to be broadly 
distributed: 



PATs{Ei,E2,Er,DQ) = 



p* 



{ei + ei + eI)Dq 



(6) 




(b) 





V(x,y) 











(a) 



O Si 



(c) 



Ba, Al 



FIG. 4: Two-dimensional representation of the plausible 
source of magnetic-field sensitive (anomalous) tunneling sys- 
tems in (e.g.) the AlBaSiO glass. The tight vitreous-Si02 
structure (a) is broken up by the Al- and large Ba-atoms (b) , 
thus leaving many metal ions free to move in a rim-minima 
(soft) tunneling potential, with n > 3 (c). The unbroken 
Si-O-Si bond dynamics, if any, is of the usual 2LS-type. 



PATsiEi, E2, E3; Dq), multiplied by the concentra- 
tion xats of these anomalous (multiwelled, and now 
near-degenerate) tunneling systems (ATS), is the proba- 
bility of finding one such ATS per unit volume. In the 
following, Xats will be absorbed in the parameter P* . 
This choice for Pats has provided a good description 
of the experimental data for the multisilicates in a 
magnetic field [l^, when in the Hamiltonian ([3]) (or 
equivalently ([S])) Dq at position is replaced with 

Dqc^'^'-' being the appropriate Peierls phase). As 
was shown in Ref. the spectrum of this B > 

modified 3LS Hamiltonian (|3]) is formally given by (using 
Viete's formula for the cubic equation's solutions): 




(7) 



-3/2 



cos 9 



which has the advantage of making use of a di- 
mensionless material-dependent parameter P*: 



(with k—0, 1, 2 and 9k = 0,+|7r,— |7r distinguish- 
ing the three lowest eigenstates) and for a choice of 
Ei,E2,E^ and Dq > ^J El + E"^ + E"^ (near-degenerate 
limit) this is shown in Fig. [5] One can see that for 
very small </> (Aharonov-Bohm phase, proportional to the 
magnetic field B) the spectrum consists of an isolated 
near-degenerate doublet which is well separated from the 
higher excited states. We shall exploit the </> = limit 
of this description for an explanation of the composition- 
dependent experiments. 
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Energy Levels, E1+E2+E3=0 




E 
b 



1 1.5 2 

Aharonov-Bohm Phase 
FIG. 5: Variation with the magnetic Aharonov-Bohm phase 
of the energy spectrum (units Do = 1) for a choice of 
El, E2, E3 with D/Do = 0.01. In this work, we are interested 
in the (p = hmit of this spectrum, which can be treated, at 
low-T, as that of an effective 2LS. 



It should be stressed at this point that in the ab- 
sence of a magnetic field, like in this work, one could 
make use of a 2LS minimal model for the description 
of the ATS, H^^\Ei,E2;Dq), and with the distribu- 
tion P{Ei,E2] Do) = P* / Da^/ El + El ensuing from the 
proximity of RER or incipient micro-crystallites. It was 
shown in Ref. |l6| that, at least for the heat capacity, this 
leads to the same physics as obtained from the 3LS mul- 
tiwelled model. There is no harm in using, for the ATS 
nesting in the incipient crystalline regions, a more realis- 
tic minimal generic multiwelled model like the above 3LS 
(3) 

Hamiltonian Hq which better approximates the physi- 
cal reality of the energy landscape. Moreover, the model 
for the composition-dependent effects remains the very 
same used for the magnetic effects and many results al- 
ready obtained for that theory can be exploited by setting 
simply B = 0. We remark, also, that a distribution of the 
type ^ for the energy asymmetry was already proposed 
for the explanation of low-T experiments with mesoscopic 
Au and Ag wires [ssl], where TS (of standard 2LS type) 
were advocated and where the poly-crystallinity of met- 
als must be accounted for. 

In summary, we have fully justified the extended TM 
which we have used in Ref. |16j and which we exploit 
also in this paper. The realistic glass is recognized to 
have a structure resembling that of chocolate [37| and as 
is pictured in the cartoon in Fig. |6l a homogeneously- 
disordered networked solid in which (at low-T in the 
glass) only standard 2LS are present with their own con- 
centration X2LS E^nd in which incipient crystallites are 



embedded. (For chocolate, these would be sugar crys- 
tals). In the proximity or within these crystallites are 
nested the ATS, with their own concentration xats 
the solid and with their own quantum mechanics and 
statistics defined by the minimal generic model repre- 
sented by Eqs. ^ and (O. This is by no means an 
ad-hoc model, since the very same model would describe 
TS in all types of real metallic and non-metallic glasses 
and quantitatively explains all of the low-T experiments 
in non-metallic glasses tackled so far. 




FIG. 6: A 2D cartoon of the chocolate-like, ceramic-glass 
structure of a real glass, in which partial devitrification has 
occurred, with the location of its low-T, two-species TS. In 
the randomly-networked bulk of the material sit the STM- 
2LS, with their own concentration X2ls, whilst within and in 
the proximity of the incipient crystaUites nest the ATS, with 
their own bulk concentration xats, each being described by 
Eqs. ([3]) and ((6]). We expect xats < X2ls and that xats 
in the best glasses. 



III. PREDICTIONS FOR THE DIELECTRIC 
CONSTANT 

The 2LS-STM has been successful in the semi- 
quantitative explanation of a variety of interesting ther- 
mal, dielectric and acoustic anomalies of structural 
glasses at temperatures T < 1 K 12, [l3|, the physics 



< 1 K 

of cold glasses being important not only for its universal- 
ities, but also because of its link with the physics of the 
glass transition (see, e.g., [1^). Beside the linear in T 
behavior of the heat capacity Cp, it is believed that the 
linear in ± In T behavior of the real-part of the frequency- 
dependent dielectric constant e'(T, w) represents a cogent 
characterization of the glassy state at low temperatures. 
We begin by deriving this behavior and putting it to the 
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test on data for e' for pure amorphous silica (no measur- 
able ATS effects, that is). 

In the presence of an applied electric field F we must 
add the dipole energy — F • po to the parameter in 

the expression ([IJ for the low-energy Hamiltonian Hq 

We can express the permittivity as e = 



F=0 



where f{F) — —-j^^ln Z{F) represents the free energy 
per unit volume. The statistical average implies also an 
integration over the two parameters of the 2LS, A and 
Ao, according to the distribution given by Eq. ([2]). We 
can write the partition function in terms of the energy 
levels £"1,2: Z = e-^^l''^'^ + ^-SiI^bT ^ 

Fig. [7] (inset) shows the behavior of the T-dependent 
part of e'(r,c^), Ae'/e' = [e'(r) ~ e'(ro)]/e'(ro), (where 
To(w) is a characteristic minimum) for pure vitreous Si02 
(Spectrosil). It can be seen that linear regimes in — InT 
for T < Tq and +lnr for T > Tq are observed, and 
roughly with slopes — -~2S and S+ = +S > 0, or in 
a -2:1 ratio. According to the 2LS-STM, in fact, we have 
the expressions 13, It 



Ae' 


Ae' 


H 


Ae' 


e' 


2LS ^ 


2 RES 




Ae' 








2RES 3eoer Jz^ 


z 


Ae' 
e' 


. PS, 

2REL 3eoer 




^/ 

J Z 


" L 







(8) 



2R.EL 



Ao,i 



2kBTz 



■ cosh (z) 



tanh ; 



1 

1 + Cj2t2 ' 



where we neglect (for low uj) the frequency-dependence 

in the RES part, where Zmm^max = ^Ormn,max/2kBT 

and where r is the phenomenological 2LS relaxation time 
given by (with E = 2kBTz) 13] 



i;AQ/7tanh 



E 



2kBT 



(9) 



In these expressions, Aomm and Aomoa; are Aq's phe- 
nomenological bounds, 7 is an elastic material parame- 
ter of the sohd and t"^^^ — _E'^/7tanh ( 



max 



^2kBT J 

i?A§^j„/7tanh ^2fc^) • P (containing the 2LS volume 
concentration, X2ls) is the probability per unit volume 
and energy that a 2LS occurs in the solid (it appears in 
Eq. ^) and p§ is the average square 2LS electric dipole 
moment. Moreover, the strategy of dielectric relaxation 
theory has been adopted, whereby the full complex di- 
electric constant e{T,Lu) has been written as, for lot <^ 1 



^RES 



(T) 



^REL 



1 



(10) 



Si02 (Spectrosil) 

2LS TM 




0.001 0.01 0.1 1 

Temperature (K) 

FIG. 7: Dielectric signature of pure a-Si02 (inset) and Al- 
BaSiO (main) glasses. Si02 data 41], fitted with Eq. ((HI, 
display a -2:1 2LS TM behavior. " AlBaSiO data 43] dis- 
play rather a -1:1 behavior, yet could be fitted with Eq. ^ 
(dashed line) [i^ with a large Aomin =12.2 mK 2LS tunnel- 
ing parameter. We have fitted all data with a more realistic 
Aomin =3.9 mK and best fit parameters from Table |T] using 
Eqs. (|S| and (fT^ (driving frequency lu =1 kHz). 



the subfixes RES and REL referring to the zero 
relaxation-time resonant and, respectively, relaxational 
contributions to the linear response e' at zero frequency. 

Presently, from expressions (|H]) we deduce that: 1) The 
so-called resonant (RES) contribution has the leading be- 
havior 



Ae; 
e' 



COCr \ 



2kBT 

Aomax 



2RES 







if T < 
if T > 



2kE 



"'(11) 



2kE 



2) the relaxational (REL) contribution has, instead, the 
leading behavior 



Ae' 







2REL 



3 CQCr \ Ao„ 



if LOTrmn > 1 
if UJTmin < 1- 



(12) 



Thus, the sum of the two contributions has a V- 
shaped form, in a semilogarithmic plot, with the mini- 
mum occurring at a Tq roughly given by the condition 
wr„.„(fc_BT) ~ 1, or kBTo{uj) ~ (^70;)^/^. eoe^ is here 
the bulk of the solid's dielectric constant and we see that 
a -2:1 characteristic behavior is justified by the STM with 
the T > To slope given by 5* = Ppl/3eoer- 

This behavior is observed in pure a-Si02 [ilj (Fig. [7] 
(inset), with the fitting parameters of TableU x = 0, from 
our own best fit to Eq. ([5])). However, in most multicom- 
ponent glasses one more often observes a V-shaped curve 
with a (roughly) -1:1 slope ratio. Fig. [7](main) shows this 
phenomenon for the multisilicate AlBaSiO glass (in fact, 
a MAS-type ceramic-glass), which has been extensively 
investigated in recent times due to its unexpected mag- 
netic field response 

0, i, H Ell ■ Also, Fig. M shows the 
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remarkable behavior of the dielectric constant vs. T for 
the glasses of composition (8102)1-2; (KgO)^; containing a 
molar concentration x of alkali oxide ;20i . It is seen that 
a 5'_/S'_|_ slope ratio of roughly -1:1 is observed, with the 
slope definitely changing with x (and faster for T > To)- 
These data from the Anderson group [2^, thus far un- 
explained by the 2LS-STM, call for an extension of the 
accepted STM and we show below that a simple expla- 
nation can be given in terms of the very same ATS that 
have been justified in Section II and advocated by one of 
us in order to explain the magnetic resp onse of AlBaSiO 
and other multicomponent glasses [16[. In view of the 
interest for these materials in low-T metrology, and on 
fundamental grounds, such explanation appears overdue 
to us. Moreover, "additional" TS (beside the standard 
2LS) of the type here advocated were already called for 
in [20J and other theoretical papers (42j . 



anomaly from the advocated ATS: 
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Here we have again neglected, for low-cj, the frequency- 
dependence in the RES part, we have put y = D/Dmm 
and TArnax is the largest phenomenological ATS relax- 
ation time given by 45 1 



"^Amax 



DVrtanh 
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2kBT 



(15) 
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FIG. 8: Dielectric signature of mixed (Si02)i-a;(K20)a: 
glasses as function of T and x [20I ]. Fitting parameters from 
Table [T] using Eq. ((8| and p4l) from our theory (driving fre- 
quency tj =10 kHz). 



For the multiwelled (3LS, in practice) Hamiltonian ^ 
we have — 3 low-lying energy levels, with < E\ <^ 
£2- 



In the E, Q a.nd D = y/Ef + + < Do 
limits (due to the chosen near-degenerate distribution, 
Eq. ([6|)) we can approximate the riu; = 3 - eigenstate 
system through an effective 2LS (though sensitive to 
all three well-asymmetries and their distribution) having 
gap AS =£1- £0: 



limA£: 



E! + Ei 



El = D. 



(13) 



(we have also exploited the condition Ex + E2 + E-^ = 
0). Using the theory of [3l to work out the 3LS- 
contributions to and e'^^;^, we arrive at the fol- 

lowing expressions for the contribution to the dielectric 



Moreover Dmin is the lowest energy gap of the multilevel 
ATS, r is another appropriate elastic constant and P* is 
the (slightly renormalised) probability per unit volume 
(after inclusion of xats) that an ATS occurs within the 
NM-pockets and channels, with p\ the average square 
ATS dipolc moment. P* and P* are so related: 



P* = P* In 



Omaa; 



Qmin 



(16) 



Domin and Domax being Do's lower and upper bounds, 
respectively. This description is intimately linked to the 
chosen distribution function, Eq. for these ATS 

which is favoring near-degenerate energy gaps D bound 
from below by Dmin- In turn, this produces an overall 
density of states given by ([11], for B — Q): 



g{E) = g2LS+9ATs{E) 



2ttP* 

2P+^^e{E-Drmn) (17) 



and that is now roughly of the form advocated by Yu 
and Legget [2l| and by some other preceeding Authors 
(e.g. [46|) to explain anomalies not accounted for by the 
standard 2LS TM. 0{x) is the step function. 

Manipulation of the expressions in ((T?| shows that: 1) 
The RES contribution from the ATS has the leading be- 
havior (note that for T < Dmin/2kB, e'\ARES is roughly 
a constant) 



Ae' 







if T < ^ 



ARES y Qeoe^kBT ) " J > 



2kB 
2kB 



(18) 



2) the REL contribution is, instead, characterised by the 
leading form 



Ae' 
e' 
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In 
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if LOTA' 
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max > 1 

max < 1. 



(19) 
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Thus, the V-shaped semilogarithmic curve is somewhat 
lost. However, adding the 2LS (Eq. and ATS (Eq. 
Pl)) ) contributions together one does recover a rounded 
semilog V-shape with a slope S'_ ~ —25* basically un- 
changed for T < Tq and an augmented slope 5+ — 
S + Sats for T > Tf) with Sats = 7TTP*pl/6eoerkBT 
that for T < Dmin/kB may approach 25* and thus (qual- 
itatively) explain a -1:1 slope ratio. 




Temperature (K) 

FIG. 9: Dielectric signature of mixed (Si02)i-a;(K20)a; 
glasses as function of T and cj for x —0.2 20]. Fitting param- 
eters from Table H] using Eq. ((8]) and (|14|) from our theory. 
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TABLE L Extracted parameters for the glasses; K-Si stands 
for the (Si02)i-x(K20)a; glasses. In all of the best fits we 
have employed the values Aomin =3.9 mK and Aomax =10 K 
extracted from fitting the pure Si02 data of Fig. [7] (inset). 

We have fitted the full expressions ([5]) and ((Ti)) to 
the data for AlBaSiO in Fig. [7] (main) and to the x- 
dependent data for (Si02)i-2:(K2 0)a; in Figs. [5]and|ni 
obtaining in all cases very good agreement between the- 
ory and experiments [20j|. Fig. IH] shows the fit of our the- 
ory to the frequency-dependent data for x = 0.2. In all 
of these best fits we have kept the value of Aomm — 3.9 
mK fixed, as obtained from our pure Si02 fit, and the 
value of Dmin also independent of x and uj. The idea is 
that these parameters are rather local ones and should 
not be influenced by NF/NM dilution. Table U gathers 
the values of all the (2LS and ATS) parameters used for 



our best fits and Fig. [10] shows the dependence of the 
prefactors (containing X2ls in P and xats in P*) with 
X. It can be seen that, as expected, the ATS prefactor 
AATfi = 'kP*p\ I (.QtrDrain scalcs linearly with x, an excel- 



lent confirmation that the "additional" TS of 2^ 42 1 are 
those ATS, by us proposed and modelled as 3LS, forming 
near and inside the microcrystallites that may nucleate 
within the NM-pockets and channels. It can be seen, in- 
stead, that the 2LS prefactor An^s — PVqI^q^t of our 
fits also increases, though less rapidly, with increasing 
X (a decrease like \ — x would be expected). We pro- 
pose (adopting a NF-, NM-cluster percolation picture) 
that new, "dilution-induced" 2LS form with alkali mix- 
ing near the NF/NM interfaces of the NF percolating 
cluster(s) as x is increased from 0. This reasoning leads 
to the expression A2LS — ^fcM/*:(l — x)+ AgurfX-^ for the 
2LS prefactor, with Ahuik, Asurf and / fitting parame- 
ters. Our best fit leads to the value / = 0.81, in fair 
agreement with the euristic expression 
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FIG. 10: The 2LS and ATS dimensionless prefactor param- 
eters (xl05) for all glasses (from Table |l]) as a function of x. 
Our data fit well with our theoretical expectations (full lines). 



f=l-iD-Ds)u 



(20) 



(where D is the fractal dimension of the percolating clus- 
ter, Dg with Ds < D is that of its "bridging" surface (not 
necessarily the hull) and is the connectedness length's 
exponent) that one would deduce from elementary frac- 
tal or percolation theory (see, e.g. [i3|)- Ds is the fractal 
dimension of that part of the NM random-cluster's sur- 
face where formation of TS takes place and we expect 
2 < Ds < D. It is indeed reasonable to expect new TS 
to be forming at these NM/NF random interfaces, for 
these are surfaces of chemical discontinuity in the mate- 
rial. The above expression is derived as follows. Imagine 
(as is shown in the cartoons in Fig. [TT|) the NM-clusters 
percolating through the NF-bulk with a site concentra- 
tion X, so that their volume scales like V ~ where 
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x'^ is their typical linear size. The number of 2LS on 



the surface of these clusters will scale like N, 



2LS 



Ds 



rf 



and so their density like N^ls/V - xx'-^^-^f'^" 
with the given expression, Eq. (|20p . for /. If wc consider 
clusters of 2D percolation and assume Dg — Dh = 7/4 
(the fractal dimension of the hull of the spanning clus- 
ter), then with D = 91/48 and = 4/3 [43| we would 
get / = 29/36 = 0.8055. More realistically, on the as- 
sumption of percolating 3D NM-clusters in the mixed 
glasses, we can make use of the values [13, El] D ~ 2.52, 
Dg = Dh = 2.14 and v ~ 0.88 to arrive at the value 
/ = 0.67 using Eq. ([20]) It is however not at ah 

clear where, at the NM/NF fractal interfaces, the new 
2LS will form (i.e. what the exact definition of Dg ought 
to be: hull surface sites, screening sites, dead-end sites 
...). If all of the hull sites are involved, then for 3D x ~ Xc 
percolation Dg = D and one then expects / = 1. Thus, 
this new phenomenology opens a tantalizing new investi- 
gation avenue for research on the applications of fractal 
theory to low-T physics. At the same time, knowledge 
of which type of NM/NF fractal interface sites are in- 
volved in the TS-formation would greatly improve our 
understanding about the microscopic nature of the TS 
(see also Ref. 23]). 




FIG. 11: (color online) A cartoon of the fractal (presum- 
ably percolating) geometry of the NM-pockets and channels 
(green), these NM-clusters growing with increasing x. 



with fit parameters from Table |n|, but it fails for the 
multicomponent glasses, like AlBaSiO, BK7, Duran (see 
e.g. [3| and references therein) and for the mixed glasses 
(Si02)i-.(K20), for :e> Q. 

Typically, the heat capacity's experimental data for 
the multicomponent glasses in zero field denote a kind of 
"shoulder" at intermediate-low temperatures. This sug- 
gests a density of states, for at least some of the indepen- 
dent TS in the glass, of the form g{E) cx l/£', in contrast 
to the standard 2LS-TM prediction, g{E) ~ const., which 
ensues from the standard TM distribution of parameters. 
Indeed, this 1/E contribution to the DOS was the very 
first observation that has led to the hypothesis of the 
ATS formulated in 

To find out the precise expression for the heat capac- 
ity due to the ATS we make use of the 3LS formulation 
for the ATS described in 16| and in more detail in Sec- 
tion II. The heat capacity is determined from the second 
derivative of the free energy with respect to temperature: 



ATS 



(T) = -T 



d'' FATs{T) 

QT2 



(22) 



where Fats {T) is the free energy of the ATS given by, if 
we neglect the third, highest energy level in the spectrum 
of Hamiltonian (|3]) (effective 2LS approximation): 



Fats{T) 



knTln e ^ 



-ksTln ( 2cosh f — ^ 



(23) 



with E = £i — Sq. The heat capacity is then obtained 
by averaging over the parameter distribution, or, equiv- 
alently, by a convolution with the DOS: 



IV. PREDICTIONS FOR THE HEAT CAPACITY 

We now come to the explanation of the, also rather 
anomalous, heat-capacity data for the mixed glasses 
(Si02)i-a;(K20), reported in [20] as a function of T and 
for different x. The heat capacity's low-temperature de- 
pendence in zero magnetic field is, for pure glasses, usu- 
ally given by the following expression: 



{T) — BphT^ + B2LsT. 



(21) 



The first term accounts for the Debye-type contribution 
from the acoustic phonons and dominates above 1 K, 
the second term is usually attributed to the low-energy 
excitations specific of all vitreous solids - the tunneling 
2LS. Bph and i?2L5 are material-dependent constants. 
This expression describes well the experimental data for 
pure silica glass at zero field (Fig. [T^l black circles: x=0 



/>oo 

C^^^'iT) = ks / dEgATs[E) 
Jq 



E 



2kBT 



cosh 



E 



2kBT J ' 
(24) 

where density of states gATs{E) has the following form 



gATs{E)^ JdD J dDoP{D,Do)5{E~D) 







if < Djnin ; 



(25) 



and Drain is the lower cutoff. The final expression for the 
ATS heat capacity results in fl6(: 



Cf'iT) = Bats 



In 2 cosh 



\2kBT 



'2kBT \2kBTj 



(26) 



13 



where the prefactor for the ATS is Bats = 
2TTP*kBnATSP{x), P* as in Section III, uats being the 
ATS mass concentration, p{x) the glass' mass density. 
Of course, xats = ^ATSPix)- For ksT > Dmin this is 
indeed roughly a constant and gives the observed "shoul- 
der" in Cp{T) when the contribution BphT^ (from virtual 
phonons) as well as the STM linear term B2lsT are taken 
into account. 

Both prefactors, for the 2LS and ATS contributions, 
are dependent on the molar concentration x of alkali- 
oxide, just as we found in Section III for the prefactors 
of the dielectric constant: B2LS — Bi,uik{^ — x)+BsurfX-^ , 
Bats — Bx. Also Bph requires to be re-evaluated. 
With increasing K2O molar concentration x for the 
(Si02)i-a;(K20)a; glass, the number of phonons from the 
NM-component (K2O in this case) increase linearly with 
the concentration x, and for the NF-component (Si02) 
it should also decrease linearly, like {1 — x). Just as 
we assumed in the previous Section III, there are frac- 
tal/percolation effects between the NM- and NF-clusters, 
which makes room for some percolation clusters' inter- 
faces where the phonons also might contribute some- 
how with a term proportional to CphX^ {Cph being an 
a;-independent constant). 

For these glasses, moreover, a non-negligible concen- 
tration of Fe'^"'" (or, according to coloring, Fe^"*") impuri- 
ties is reported, a side effect of the industrial production 
process. Estimates give 102 ppm for AlBaSiO and 126 
ppm for Duran, 6 ppm for BK7, 100 ppm for Pyrex 7740 
and 12 ppm for Pyrex 9700 (see, e.g., the discussion in 
Rcf T|)- All glasses may, indeed, contain some [Fe04]'^ 
impurity-substitution F-centers (in the glass, similar to 
a liquid, in concentrations however much, much lower 
than the nominal Fe bulk concentrations [341 )• The Fe'^"'" 
cation and the anion, on which the hole is localized 
(forming the 0~ species, that is the -I- hole subsys- 
tem), form a bound small polaron. In this configuration 
the Fe'^"'" cation is subject to a crystal field with an ap- 
proximate C3 symmetry axis along the Fe'^"'' - 0~ direc- 
tion. This axis plays a quantization role for the Fe'^"'" elec- 
tronic spin. The hole is assumed to be tunneling between 
two neighboring oxygen ions, switching the quantization 
axis between two directions and therefore entangling its 
spin states. This is likely to give some tiny contribu- 
tion to the heat capacity and we should, therefore, also 
take it here into account f50| . The spin Hamiltonian of 
the [Fe04]° F-center is Hg-s = V^s^Sz, where Vz is the 
principal value of the dipole interaction matrix, Sz and 
Sz are the spin operators of the hole and of the Fe'^^ 
ion, respectively. In the absence of a magnetic field there 
are only two low- lying energy levels: £^1^2 = ±| |^z|- The 
unknown distribution function G{Vz) must approach zero 
when its argument approaches either zero or infinity and 
have a maximum at a definite argument value Vq. The 
simplest one-parameter function displaying such proper- 
ties is a Poisson distribution: 



G(14) = -^expf— ^j, Ke(-oo;Q], K) < 0. 

(27) 

The contribution from the [Fe04]'' ensemble to the heat 
capacity is, as usual: 

C^,3+(T) = -T^^g|^, (28) 

where Fpe3+ (T) is the free energy of the [Fe04]° ensem- 
ble, that one evaluates as: 

= Fp,,+ = -keT In (^2cosh ) , (29) 

(here E — j\Vz\)- Using said distribution function for 
G{Vz), Eq. (P7)) . as well as the expression for Cp^i+ {T) 
from Fpf,3+ (T), one can obtain an expression for the heat 
capacity from the trace [Fe04]° centres in the glass, and 
which should be added to the total heat capacity Cp: 

Cp"'* {T) ^ p{x)njkB X 

X HdYzf -^Vcosh-2 ( G{Vz) 

^pix)n,kBfdVz'^^e^-'^hosl.-^ (H^) ' 

(30) 

where Uj = Xj/p{x) is the mass concentration of the 
tiny amount of Fe^"*" ions (a very small fraction of the 
total bulk Fe-concentration) substituting the Si''"'" in the 
network. 

Hence, the total heat capacity will be the sum of all 
these contributions (Eqs. jtl]), (jUl) and ^) jSll]: 

Cp{T) = BphT^ + B2LsT + Cf^{T) + C^^"^ (T). (31) 

Making use of expression ([SI]) , we have fitted the 
experimental data for the heat capacity of the 
(Si02)i-a;(K20)2, glasses from [2^. In order to fit the 
pure a-Si02 data we use only formula (f2T|) . that fits the 
pure silica's data well within the 2LS-STM. 

The heat capacity Cp{T,x) data ^o[ for the 
(Si02)i-a;(K20)2: glasscs were obtained using a signal- 
averaging technique and for these samples the data are 
presented in Fig. [121 As one can see, the heat capacity 
for the (Si02)i-a;(K20)a; glasses at low temperatures is 
larger than that for pure silica glass, as is typical for the 
multicomponent glasses, already with the smallest 5% 
concentration of K2O. The heat capacity decreases and 



14 



then again increases with increasing molar concentration 
X of K2O. The additional heat capacity arises from the 
addition of ATS in the K2O NM-clusters and also from 
the presence of Fe"^"*" impurities, contained in small (and 
unknown) concentrations, but contributing to the low- 
and middle-range of the temperature dependence. 




10"' 10° 
Temperature (K) 

FIG. 12: The temperature dependence of the heat capacity 
for a-Si02 (black circles) and for the (Si02)i-i(K20)2, glasses 
[20t . The full lines are our theoretical curves, as generated by 
Eq.dSIl). 



glass 


X 


Bph X 10* 


B2LS X 10** 


Bats x 10* 




type 


mol 


Jg-^K-^ 


Jg-^K-^ 




ppm 


SiOa 





245.55 


70.65 






K-Si 


0.05 


260.92 


155.23 


22.77 


29.86 


K-Si 


0.08 


266.36 


196.11 


36.44 


18.15 


K-Si 


0.10 


269.46 


221.62 


45.55 


10.54 


K-Si 


0.20 


281.42 


337.19 


91.11 


3.00 



TABLE II: Extracted parameters for fits to the heat capacity 
data for Si02 and (Si02)i-2;(K20)a: glasses, with Dmin =0.87 
K and Vo =-0.42 K as fixed. 

Both prefactors, for 2LS and ATS, are indeed de- 
pendent on the molar concentration x from our data 
analysis, and in the same way as we did in Section III 
we have fitted the extracted prefactors with the forms: 
B2LS - Bbuik{i - x) + BsurfX^, Bats - Bx {B being 
some constant). These dependencies are show in Fig. [T31 
Also Bph is found to change by increasing the concentra- 
tion X of the good crystal-former, K2O, and in the way 
we anticipated. 

With increasing concentration x, for the 
(Si02)i-x(K20)2^ glass the number of phonons from the 
NM-component (K2O) increases linearly with the con- 
centration, and for the NF-component (Si02) it should 
be decreasing linearly like {1 — x). As we reasoned for the 
dielectric constant, there are percolation mixing effects 



f 

^ph=^Si02<^-'''"'^K20'<*Cp|^x 
B„;„,=245.5S, B^,n=0, C„h=312.95 


■ 




■ Phonons 




• 2LS 




* ATS 


'2LS=Bb(1-'')+Bs'<' 




Bg=1 047.29 










B=455.55 
■ 



0.00 0,05 0.10 0.16 0.20 

Concentration x 

FIG. 13: The 2LS and ATS prefactor parameters (x 10*) for 
all glasses (from Table HI)) as a function of x. The experimental 
data fit well with our theoretical expectations with / = 0.81 
(full lines). 

between the NM- and the NF-systems, which create 
percolation clusters and their NF/NM interfaces where 
phonons also might be populated in a way proportional 
to CphX-^ . As it turns out, the very same value / = 0.81 
can be extracted from all our fits, just as was done in 
Section III for the dielectric constant data. 



V. SUMMARY AND CONCLUSIONS 

Thus, we have shown that there is direct evidence 
in zero magnetic field already for the existence of mul- 
tiwelled ATS (modelled as tunneling 3LS) and with 
the new distribution function advocated to explain the 
magnetic-field effects in the multicomponent glasses (see 
[lit). The relevance of near-degenerate multiwelled TS 
in glasses is a new and unexpected finding in this field of 
research. Our work predicts, in particular, that the mag- 
netic response of the mixed alkali-silicate glasses should 
be important and scale like the molar alkali concentra- 
tion X. At the same time the -1:1 slope-ratio problem of 
the standard TM in comparison with experimental data 
for e'(T) has been given a simple explanation in terms of 
our two-species tunneling model. The main result of this 
work is that the concentration xats (absorbed in P* and 
thus in the Aats- and B^rs-prefactors) of ATS indeed 
scales linearly with x for both e'(T, x) and Cp{T, x). This 
is supported by our analysis of the experimental data of 
Ref. [20[ very well indeed. This analysis is, in our view, 
strong evidence that the ATS are nesting in the NM- 
clusters of the good crystal-formers. 

Using the results of this analysis (and for AlBaSiO 
the results of the experimental-data analysis in a mag- 
netic field we can estimate the value of the dipole 
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ipj. For Al- 



moment associated with the ATS, Peff 

BaSiO, using the value of P* extracted from Cp [16j and 
that of Aats given in Table H] we extract Pe//=0.41 D. 
For (Si02)i-K(K20)2;, we notice from the definitions in 
Section III that the ratio of the dielectric and heat ca- 
pacity prefactors, 



A 



ATS 



p{x) 



Bats 2eQerkBDrmnPi'^ 



(32) 



is almost independent of the K2O concentration x. From 
our extracted values in Tables U and and the mea- 
sured values of p{x) [l^ we estimate pe// — 0.045 D for 
the mixed glasses, independently of a;! Considering the 
elementary atomic electric-dipole's value is eap = 2.54 
D, these small values of Peff for the ATS confirm that 
their physics must come from the coherent (or correlated) 
tunneling of small ionic clusters (the very same origin 
for the large values of Dmin and for Domin,max, see the 
Appendix). Indeed, a cluster of N coherently-tunneling 



Ei=lPi 



that can 



particles has a dipole moment Pef / 

become much smaller than eao (the order of magnitude 
of each jp^j in the sum) as N grows large. The fact, that 
we extract values oipeff much smaller than eao, confirms 
the picture of a correlated tunneling cluster in the B = 
case already. 

It is noteworthy that several papers from the Ander- 
son group have proved that the addition of any NM- 
species in a networking pure glass causes significant (and 
thus far unexplained) deviations from the predictions of 
the 2LS-STM [lO, [Hj. We have explained the origin of 
these deviations, for Cp(T, x) as well as for e(T, x). How- 
ever, experiments do show that the thermal conductivity 
k(T,x) cx remains (below 1 K) remarkably univer- 
sal and composition-independent [20[. This is connected 
with the super-universality of the internal friction coeffi- 
cient, Q~^, in the cold glasses; these and other remark- 
able findings will be addressed elsewhere within the con- 
text of our approach. 

In summary, we have shown that there is direct ev- 
idence in zero magnetic field already for the multi- 
welled ATS advocated to explain the magnetic field 
effects in the multicomponent glasses. Similar x- 
dependent phenomena are to be expected for the low- 
T anomalies of the MAS-type ceramic-glass of composi- 
tion (Si02)i-a;(MgO)a;, which should also respond to the 
magnetic field [52| (just like the mixed alkali-silicates of 
this work should). One may remark, at this point, that 
any extension of the 2LS-STM enlarging the adjustable- 
parameter space is bound to improve agreement with the 
experimental data. In this paper we have shown that 
it was not just a matter of quantitative agreement, but 
qualitative as well. Whilst agreeing that the TM remains 
unsatisfactory, we stress that it is the only approach we 
know of that is versatile enough to allow for an interesting 



explanation of rather puzzling phenomena at low-T in the 
real glasses. Furthermore, our two-species, multilevel TS 
model has been able to consistently explain a good num- 
ber of different experimental data 3, 3] ■ It cannot be a 
mere coincidence that the same phenomenological model, 
with rather similar material parameters in different ex- 
periments, is capable of explaining so much new physics. 
Far from being an ad-hoc model, our approach reveals 
the intimate microscopic structure of the real glasses, 
which cannot be considered as being homogeneously- 
disordered anymore, and this must have some important 
consequences also for a better understanding of the mech- 
anisms underlying the glass transition. 

As for the possibility of estimating the size and den- 
sity of the incipient crystals in glasses from our theory, 
we remark that the simplified geometric-averaging pro- 
cedure adopted for the physics of the ATS so far [3] 
does not allow anything more than an estimate of the P* 
parameter (- 1.97 x lO" cm^^ fo^. AlBaSiO this 
being in fact the value of xatsP* ■: P* being the unknown 
dimensionless parameter of the ATS distribution in Eq. 
dni). However, the geometric-averaging procedure should 
be performed in two stages (within the incipient micro- 
crystals first and then within the glassy matrix in which 
the crystallites are embedded) at the price of making 
the theory considerably more complicated. When this is 
done, with a more efficient and complete theoretical for- 
mulation, then information on the size distribution of the 
incipient crystallites could be gained from further low-T 
experiments in magnetic fields and at different controlled 
compositions. 
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APPENDIX 

We first show how the spectrum of the 4LS Hamilto- 
nian, Eq. ([4]), is similar to that of the 3LS Hamiltonian, 
Eq. ^ in the near-degenerate limit. We rewrite the h'^^ 
4LS Hamiltonian in the presence of a magnetic field, cou- 
pled orbitally to the charged tunneling particle: 
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^(4) _ 



( El Die'-^/^ £,^(,^'1^/2 D^f.-i'P/i ^^ 
Die-'-'I'l'^ E2 Die*'*/'* Dae*'*/'* 
£)2e-»"*/4 Die""*/"* Die*'*/" 
V Die*'*/'* Dae-*'*/^ Die~*'*/* D4 / 

(33) 

where = is imposed, Di and D2 are the n.n. 

and n.n.n. hopping energies, respectively, and where 



2tt 



£(B) 



$(B) = So-B, $0 



he 



(34) 



is the Aharonov-Bohm phase resulting from the magnetic 
flux $(B) threading the square-loop (having area 5*0) 
closed trajectory of the particle. The above Hamilto- 
nian should in fact be symmetrised over its permutations, 
since the sign of the n.n.n. Peierls phase is ambiguous (in 
practice, one replaces Dae^*'*/^ with D2COs(0/2) in the 
appropriate matrix entries). The eigenvalues equation 
giving the energy levels is then: 



£^+£^(y^ E,Ej - 4D'i - D'i{l + cos < 



£( } E^EjEk+4DlD2{l+ COS c 



(35) 



y EiEj 

i<j<k 

+E1E2E3E4 — D\{EiE2 + D2D3 + E3E4 + E4E1) 
-^DliEiE^ + E2E4){1 + cos(j)) - 2D^D|(1 + cos(? 

+2D:f (1 - cos (/)) + ^D^ (3 + 4 cos + cos(20)) = 0. 



More instructive than numerically extracting the four ex- 
act roots fo,i,2,3 (with Sq < £1 < £2 < £3) is for us the 
physically interesting limit case in which |Di/Di| <C 1 
and ID2/D1I <^ 1 (near-degeneracy of the four- welled 
potential). The above eigenvalue equation then becomes 
much easier to study: 



£'^ + £^(J2^^^:> -'^^i 



— D\{ExE2 + -B2-E'3 + £"3^4 + E1E4) 

-H2D^(1 - cos?!)) w 



(36) 



this being the eigenvalue equation of the reduced 4LS 
Hamiltonian 



H. 



(4) 

Orerf 



/ Di Die"*/'* Die-"*/'* \ 

Die"*'*/'* D2 Die*'*/'* 
Die-*'*/" D3 Die*'*/" 
V Die*'*/'* Die"*'*/'* D4 / 

(37) 



always for |Di/Di| ^ 1, and which has the following 
solutions: 



^0,1,2,3 ^ ^ J ^ EiEj 



Dl 



± 



D,D,\2 



i<3 



Dl 



+4 



E1E2 + E2E3 + E3E4 + E4E1 



+8{cos(l)~ 1] 



1/2 



1/2 



(38) 



The perhaps surprising result is an energy spectrum 
where only the middle doublet {£1 , £2 in our notation) be- 
comes near-degenerate at weak (or zero) magnetic fields 
((/) — )■ 0). This is shown in the inset of Fig. [14] and is 
reminiscent of the situation with dimcrizcd 2LS consid- 
ered in Ref. 55| in order to account for the oscillations 
of the dielectric constant with B. Beside there being 
no evidence for a dimerization of TS in glasses (unlike 
perhaps in mixed and disordered crystals), one would 
have to explain why the ground state £0 is prohibited 
for the tunneling particle (the real energy gap being in 
fact A£ = £1 — £q). The way out can be found again 
in Sussmann's paper (36j since the n^, = 4 welled trap- 
ping potential giving rise to the same physics as our 3LS 
must in fact have tetrahedral and not square geometry. 
The tetrahedral 4LS in a magnetic field will be consid- 
ered elsewhere Here we only want to remark that 
the tetrahedral situation can be mimicked by a square 
multiwelled potential in which jDi/Da] <^ 1 and always 
in the limit case jDi/Dal ^ 1. This corresponds to the 
counter-intuitive situation in which it is easier for the 
particle to tunnel across the square to the n.n.n. site 
rather than to a n.n. site, as if the middle potential bar- 
rier had collapsed. In this limit case, Eq. (1551) becomes, 
instead: 



£'^ + £^(Y^ EiEj - 2Dl cos^ - 



-DliEiEz + E2E4) + Dl cos" I « 



(39) 



with, once more, easily found solutions 



^^ = ±^i2cos2- / , 
^2 V2I 2 ^ Dl 



EiEj 



± 



E1E3 + E2E4 . (/) 

+4 4 cos - 



1/2 



1/2 



(40) 



as exemplified in Fig. [T3] (main). We therefore ob- 
tain that the lowest-lying gap remains near-degenerate 
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for weak fields and we can conclude, therefore, that the 
lowest-lying eigenvalues display, at low-T, almost the 
same physics as in the case of a 3LS (Fig. [S]). This shows 
the important role of the frozen solid surrounding the 
tunneling "particle" (which could be, perhaps, a vacuum 
in fact) and that when nested within an incipient crys- 
tallite a magnetic-field sensitive TS is well described by 
a tunneling 3LS as the minimal generic model potential. 

Energy Levels, E1-hE2-hE3+E4=0 



potentials we expect instead: 




1 2 
Ahai'onov-Bohm Phase 



FIG. 14: (Main) Variation with tlie magnetic Aharonov- 
Bohm phase of the energy spectrum (units D2 = 1) 
for the case Di — and a choice of Ei, E2, Es, E4 with 
^yEf + E^+Ei+ EI/D2 = 0.01. This is to be compared 
with the 3LS energy spectrum. Fig. (5] (Inset) The en- 
ergy spectrum in the opposite case, D2 = 0, and a choice 
of El, E2, E3, Ei with Di about 100 times stronger. 

Next, we seek a description of a cluster of N corre- 
lated tunneling particles (atoms, ions or molecules) and 
derive the transformation rules for the tunneling param- 
eters (including also those involved in the theory for the 
magnetic effects [l^) when the cluster is replaced by a 
single "tunneling particle" as the result of the coherent 
tunneling (CT) of the particles in the cluster. 

The "tunneling particle" in question is only a fictitious 
one, as was inferred in Section II by examining local min- 
ima in the energy landscape, representing the CT of a 
cluster of N true tunneling particles (which in the real 
glasses might be the lighter species involved in the ma- 
terial: Li"*" in the disordered crystal Li:KCl, in the 
multisilicates and 11+ and/or D+ in a-glycerol) and for 
which we have to make up appropriate renormalized tun- 
neling parameters. The concept of CT in separate local 
potentials is distinct from that of the joint tunneling of 
N particles in the same local potential, for in the latter 
case the tunneling probability would be depressed expo- 
nentially: Da/n w 17(^)^ (Ao being the real particles' 
common tunneling transparency). As we shall show be- 
low, at least for moderate values of N , for CT in separate 



Do « iVAo, 



(41) 



and, for the fictitious particle's charge and fiux-threaded 
area (see [l^ for the magnetic effects): 



S 



A 



iNai 



(42) 



(go — 0(e) being the charge of the real tunneling parti- 
cles, oq Bohr's radius). In the latter relations, less obvi- 
ous is the renormalization of the fiux-threaded area 5a 
of a 3LS ATS. It is however the direct consequence of our 
multiphase model of a real glass, thought of as made up 
of regions of enhanced atomic ordering (RER) or micro- 
crystals (Figs. [21 ini) embedded in a homogeneously- 
disordered host matrix. The magnetic flux appears 
quadratically in our theory [16| , each elementary fiux 
adding up within each micro-crystallite or RER and then 
appearing, squared, multiplied by cos^ /? in the glassy ma- 
trix in a magnetic field (/? being the random angle formed 
by Sa with the magnetic field B), a factor averaging out 
to ^ in the bulk. From these considerations and from 



Eqs. pi]) and p2|). the renormalization of the composite 
phenomenological parameter Do\-^\Sa would be as fol- 
lows (if q = 2e, appropriate for the multisilicates): 



(43) 



Setting Ao=l mK, one gets a value of N ranging from 
about 25 coherent-tunneling particles in a cluster at the 
lowest temperatures [i^ , to about 600 at the higher tem- 
peratures. These estimates are somewhat speculative, 
since the real values of the elementary fiux-threaded area 
and of the elementary tunneling barrier transparency Ag 
are unknown, we are however inclined to support the 
value N « 200 that was proposed by Lubchenko and 
Wolynes [2J|. This would yield a value of Ami™ ranging 
from 80 /zK to 4 mK also for the mixed alkali-silicate 
glasses (for which I?„ii„ ~ 800 mK). The above consid- 
erations show all in all the tendency for the coherent- 
tunneling cluster size N to be also temperature depen- 
dent. 

We now come to the justification of Eq. (|4T|) . At 
low temperatures the interactions between true tunnel- 
ing particles become important and coherent-tunneling 
motion can take place. Coherent motion in the context 
of the tunneling model is a state in which all of the parti- 
cles in each local potential contribute to the overall tun- 
neling process in a correlated way. We exemplify our 
ideas in the context of the simplest 2LS situation first. 
Let us consider two interacting 2LS. Let the positions 
of the particles in the two wells be left (L) and right 
(R). The tunneling particles in the cluster interact via a 
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weak potential U which may have its origin, for example, 
from either a strain-strain interaction having the form 
U ^ A/r^ (dipole-dipole interaction) [24, 53|, where r is 
the distance between a pair of tunneling particles either 
in the L or R wells and ^ is a constant, or it could be 
due to electrostatic dipole-dipole interaction. The tun- 
neling of the particle in one 2LS from L to R (or vice 
versa) influences, via the interaction, the particle in the 
other 2LS, forcing it to jump into the free well. The hop- 
ping Hamiltonian of two interacting 2LS can be written 
as follow (with A^^ = — A^^^ = A^ and dropping factors 
of -1/2): 

H2 = ^ Aiac\^Cia + Aoi ^ c\aCla' + hc 
a—L.R a-^a' 

+ ^ A2a4aC2a + Ao2 ^ 4aC2a' + hc 

- U {c\^CiLcl^C2L + c\j^Cib.4j^C2R^ (44) 

which favors coherent LL— >RR and RR— >LL joint 
tunneling and acts on the joint states |aa') = 
\LL), \LR), \ RL), \RR). The coherent motion of the two 
real particles can now be replaced by the tunneling of a 
new, fictitious particle in its own double well. In order 
to write the renormalized Hamiltonian of two coherent- 
tunneling particles we are interested only in the matrix 
elements {LL\H2\LL), {RR\H2\RR), {RR\H2\LL) and 
{LL\H2\RR) of Hamiltonian jB]): 

{LL\H2\LL) = A1 + A2-C/ (45) 

{RR\H2\RR) = -A1-A2-C/ 

{RR\H2\LL) = {LL\H2\RR) = Aoi + A02 

(instead of the latter two, the pair {RL\H2\LR) and 
{LR\H2\RL), having the very same value Aoi + A02, 
could have served the purpose). These matrix elements 
represent the Hamiltonian of the fictitious particle, which 
corresponds to both real particles tunneling coherently 
together: 



Ai + A2 - [/ Aoi + A02 
Aoi + A02 -Ai - A2 - [/ 



(46) 



(the condition A[ + A2 = to be fixed through the addi- 
tion of an overall constant). Next, we consider the case 
of three interacting 2LS and repeat the previous consid- 
erations. The Hamiltonian of three interacting 2LS has 
the form: 



The matrix elements of a single replacing fictitious par- 
ticle that correspond to CT are obtained as follows: 



{LLL\H^\LLL) = Ai + A2 -f A3 - 3C/ (48) 
{RRR\H3\RRR) = -Ai - A2 - A3 - 3C/ 
{RRR\H3\LLL) = {LLL\H3\RRR) = Aoi + A02 + A03 

(the choice of the latter two not excluding the remain- 
ing coherent matrix elements pairs: {RRL\H3\LLR) and 
{LLR\H3\RRL), {RLR\H3\LRL) and {LRL\H3\RLR) , 
{LRR\H3\RLL) and {RLL\H3\LRR) , which are all 
equivalent). One can notice that the renormalized tun- 
neling parameter is the sum of the Api of each 2LS. The 
energy asymmetry is also the arithmetic sum of the A^ 
of each 2LS, but one must add the interaction energy 
-U multipHed by N{N - l)/2. Thus, for a coherently- 
tunneling cluster of N 2LS we find that the diagonal 
matrix element becomes, generalizing to arbitrary N: 



A = Eti A 

that corresponds to the CT-splitting for all N particles, 
becomes simply Ao = Aoi. 

Applying the previous considerations to our model for 
a number N of ATS with three wells (see Fig. |4]) we can 
write the interacting Hamiltonian in the form: 



N{N- 



—U and the off-diagonal element. 



TV I 3 

^ '^E,a4^C,a + {Dot4^C,a' + hc) 

i—1 I a— 1 a^a' 

' C^EE4^-4aC.'a (49) 



If we represent the group of N coherently-tunneling parti- 
cles as a single fictitious particle moving in a three- welled 
potential, which is characterized by its own ground state 
energies Ea and tunneling parameter Uq, we can describe 
this renormalized 3LS by the following Hamiltonian: 



= E! Eaca^ca + V D^CA^CA' + he (50) 

^-^ A=\ ^-^A^A' 

The ground states energies Ea in the wells and tunneling 
parameter Dq for the fictitious particle, in line with the 
calculations above, can be obtained through: 



Ea 
Do 



a\HN\aa ... a) 



A 



1,2,3 



(a'a' . . . a'\HN\aa . 



iaa . 



. alHNla'a 



(51) 



i=l I a=L.,R 



^ E/ E/ 4a^ia-4' a'^i' a 



i<i' a 



(and the remaining variants of the second definition line, 
all equivalent). In analogy with the 2LS considerations, 
one can sec that the renormalized tunneling parameters 
(47) D = Ef + E2 + E'^ and especially Dq can be replaced 
by the arithmetic sums of those of the bare coherently 
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tunneling particles, D « NDi (neglecting the correction 
for a sufficiently weak U and moderate values of N) and 
Dq sa NDoi, respectively. Indeed, the tunneling prob- 
abilities of weakly-correlated events should add up for 
values of N not too large. Therefore, since TV can attain 
values as large as 200 [2J] (independently of the solid's 
composition) in some models, and as corroborated by 
our reasoning in this Appendix, this leads to values of 
Di and Doi (as extracted from our theory's fitting pa- 
rameters) comparable to those characteristic of the 2LS 
TM. The large values of Dmin and especially of Domm 
and Domax , as extracted from our fits of our theory to the 
available experimental data, find therefore an interesting 
and physically plausible explanation. 




FIG. 15: A cluster of N = i weakly interacting (real) tun- 
neling particles that is being replaced with a (fictitious) single 
3LS (Fig. \M.c)) having renormalised parameters according to 
Eqs. (|4T|l and (|42|l . 
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